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Abstract 

We analyze the dynamics of a Brownian gas in contact with a heat bath in which large tem- 
perature fluctuations occur. There are two distinct time scales present, one describes the decay of 
the fluctuations in the temperature and the other one is associated with the establishment of local 
equilibrium. Although the gas has reached local equilibrium, there exist large fluctuations in an in- 
tensive parameter (temperature) which break the thermodynamic equilibrium with the heat bath. 
Thus the decay of the fluctuations in the intensive parameter is larger than the characteristic time 
for the establishment of local equilibrium. We show that the dynamics of such large and intensive 
fluctuations may be described by adopting a nonequilibrium thermodynamics approach with an ad- 
equate formulation of local equilibrium. A coarsening procedure is then used to contract the space 
of mesoscopic variables needed to describe the dynamics of the gas and the extensive character of 
the description is lost. This procedure allows us to derive an effective Maxwell-Boltzmann factor 
(EMBF) for the Brownian gas, as has been recently proposed in the literature [Ej- Furthermore, 
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we use this local equilibrium distribution and an entropy functional to derive a nonequilibrium 
probability distribution and a hydrodynamic description for the Brownian gas which contains fluc- 
tuating transport coefficients. The ensuing description is nonextensive and our analysis shows that 
the coarse-graining procedure is responsible for the nonextensivity property 

Author Keywords: Nonequilibrium systems, Fluctuations of temperature, Effective Boltzmann fac- 
tor, Tsallis statistics 

PACS numbers: 05.40.-a 
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1 Introduction 



Complex systems where long range interactions, large external gradients or large fluctuations of their 
intensive variables are present, show peculiar behaviors which have attracted a good deal of attention in 
recent years El El El Hj • The theoretical description of these systems constitutes a nowadays challenge 
and is far beyond the limit of applicability of the theory of fluctuations near equilibrium as developed 
by Onsager and Machlup [Zj and subsequently generalized by Fox and Uhlenbeck 8J. This is essentially 
due to the fact that these theories are restricted to describe fluctuations of extensive variables only. 

Several theories have been recently proposed to describe the special properties of simple and complex 
fluids in nonequilibrium stationary states El El El IU E] • However, theories describing physical situ- 
ations where fluctuations of the intensive variables, such as temperature or chemical potential among 
others, may be important are scarce, jXTJ These issues have been mostly investigated within the 
framework of the so called nonextensive statistical mechanics (NESM) [TT], ^2]; which according to 
Beck and Cohen ^J|, is only one of the possible more general statistics that can deal with the physical 
situations mentioned above. NESM proposes that for systems with a sufficiently complex dynamics, 
its necessary to introduce a more general statistics than the usual Boltzmann-Gibbs description. Such 
statistics have been discussed for a variety of systems J3| and have led to the introduction of an effective 
Boltzmann factor (EMBF) which reduces to the ordinary one in the appropriate limit. 

These concepts have been explored in more detail within the model of Brownian motion, for which 
Beck and Cohen have postulated a modification of the Boltzmann factor whose explicit form depends 
on the statistics of the intensive parameters fluctuations. They show that if it is assumed that the 
inverse temperature of ordinary statistical mechanics fluctuates with the Gamma distribution, then 
the ensuing Brownian dynamics can be correctly described by Tsallis statistics, |TH1 Elj. One of the 
main purposes of this paper is to use this same model of nonequilibrium Brownian motion to examine 
how the generalization of the Boltzmann factor may be justified in terms of well established concepts 
used in irreversible thermodynamics, such as Einstein's formula for the probability distribution and a 
Kullback-like entropy functional [13 EE1 H3 HH1- By adopting this point of view, we show that the 
modified Boltzmann's factor has its origin in the contraction of the space of mesoscopic variables that 
describe the dynamics of the system, and that this contraction is responsible for the non-extensive 
character of the new description. 

Usually, the dynamics of a Brownian gas is analyzed by taking into account only the equilibrium 
fluctuations of the extensive variables of the heat bath, and by assuming thermodynamic equilibrium 
between the Brownian gas and the bath. Since the characteristic decaying time of these fluctuations 
is short as compared with the time characterizing the establishment of local equilibrium, they become 
independent from those of the velocity of the Brownian particles which, in this case, are well described 
by the Maxwell-Boltzmann factor . However, a more general case occurs when intensive parameters 
(temperature) may fluctuate and are large enough to break its thermodynamic equilibrium with the 
heat bath. In this case, the Brownian gas can be assumed to be in local thermodynamic equilibrium, 
but not with the entire heat bath. In this situation the work performed on the gas will be different from 
that corresponding to the equilibrium case. These conditions can be met when the system is constrained 
to be far away from equilibrium and the time scale and wavelength of the intensive fluctuations are 
larger than those characterizing the establishment of local equilibrium. As a consequence, the proba- 
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bility distribution for the fluctuating velocities of the Brownian particles depends, in general, on the 
distribution of the temperature fluctuations and is no longer described by the usual Maxwell-Boltzmann 
factor. 

The paper is organized as follows. In Sec. II we define the system and obtain a probability 
distribution for the velocities of the Brownian particles that incorporates the effects of nonequilibrium 
temperature fluctuations; from it we derive the effective Maxwell-Boltzmann factor. Then in Sec. Ill 
we derive a Fokker-Planck equation describing the relaxation mesoscopic dynamics of the Brownian gas 
and we use it to derive a hydrodynamic description with fluctuating transport coefficients. A generalized 
Onsager-Machlup formalism is presented in section IV and the last section V is devoted to discuss our 
main results. 



2 Model and Basic Equations 

As mentioned above, the dynamics of a diluted gas of Brownian particles in contact with a heat bath 
in the presence of two distinct time scales, one related to the decay of the fluctuations and the other 
to the establishment of local equilibrium, can be interpreted by assuming that each Brownian particle 
is embedded in a system in contact with a heat bath at temperature Tg(i). We want to analyze the 
dynamics of the Brownian particle when nonequilibrium temperature fluctuations 5T(t) occur in the 
system. In particular, we will assume that the characteristic length and time scales of these fluctuations 
are sufficiently larger than those characterizing the establishment of local equilibrium. The relation 
between the system and heat bath temperatures is then T(t) = Xb(£) + 8T(t). 

When a fluctuation occurs, the system only performs the amount of work W on the Brownian 
particle. This modifies its kinetic energy by W/m = Sekin = — -rr, where m is the mass of the particle 
and the sign indicates that the system performs the work. Moreover, the system exerts the work p5V 
and transfers the amount of heat T5S to the bath. Assuming that the variation of the total entropy 
of the system plus bath is 5 St — SS + 5Sb, and by using the Gibbs equation, after dividing by m we 
obtain the relation [UJ [T7] 

5s T = t^t- [Se kin + 5e + pSp' 1 ] + — ^— [5e B + Pb^P'b 1 ] , (1) 
l(t) l B {t) 

where 5e and <5e# are the variations of the internal energy of the system and the bath, p and ps the 
system and bath pressures, and p and ps the corresponding specific volumes. 

Eq. (Q) can be used to obtain the minimum amount of work performed by the system to modify the 
state of the Brownian particle [14] . This situation occurs when the processes taking place between the 
system and the heat bath are reversible. Since in that case SeB = —Se and 5p B x = — <5p _1 , and given 
that 5T = T — Tb, then Eq. (JTJ) becomes 

Ss T = 7^5e kin - T^STSe + ( £ - ^ J 5p-\ (2) 
T(t) TT B \T T b ) 

which is in agreement with the Le Chatelier-Braun principle [T4"] 117] . In order to write equation (j2j) in 
a more convenient form, we expand Se(T,p) up to first order in 5T and After using the relation 
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which can still be simplified if we use again the relation T = Tg+ST in order to make the approximations 
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into Eq. (j3J), keeping terms up to second order in the 
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where Sp = p — ps- After substitution of Eq. 

fluctuations and using the expansion Sp ~ (jj^) Sp + (^)^<5T, we arrive at the following expression 
for the variation of the total entropy Sst 
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(5) 



However, since our interest in this paper is to analyze the effects of the large temperature fluctuations 
on the dynamics of the Brownian gas, we shall adopt a statistical description point of view and assume 
that the probability Pi e of occurrence of a fluctuation in the composite system may be expressed through 
Einstein's formula 14 j 

P le = P e^T, (6) 
where ks is Boltzmann's constant and Pq normalizes P\ e . Thus, by substituting Eq. (j3J into © and 



using 5e 



kin 



2 ' 



Pi P can be rewritten in the more convenient form 



P le = p e » b tbW 



(7) 



which is the probability that a fluctuation occurs. To arrive at Eq. (J7J) we have expanded 



Kj.hiST) with h(ST) 



1 - 



ST 
T B (t) 



8T 
T B (t) 



This last equation allows us to analyze how the 



distribution function accounting for the fluctuations of the velocities of the Brownian particle becomes 
modified in the presence of large temperature fluctuations ST occuring in the system. Integrating Eq. 
(0) over 5T leads to the effective Maxwell-Boltzmann factor P>i e 



B, 



le 



2k B T B (t) 



h(8T) 



d(ST), 



where we have assumed a constant volume processes for simplicity. This probability distribution is 
an effective or modified Boltzmann factor (EMBF) for the nonequilibrium Brownian system and it 
contains corrections arising from a non-Gaussian distribution. However, in the usual case where system 
and bath are in thermal equilibrium, h(5T) — > 1, and the temperature fluctuations become independent 
of velocity fluctuations, as expected [13] . 



5 



2.1 Connection with superstatistics 

Equation (jHJ) implies that when system and bath are not in thermal equilibrium, temperature fluctua- 
tions modify the local equilibrium distribution of the Brownian particle in the sense discussed in Ref. 
[T3*] . To illustrate this point more explicitly, let us rewrite (jHJ) as 



Bie(e B ) 



f(5T) e- h{ST)eB d{5T) = (e~ h{ST)eB ) f , 



(9) 



where e B = ^f^y 
fluctuations / (ST) given by 



and (..)/ denotes the average over the distribution function of the temperature 

(10) 



mh(ST) 
k B T B (t) 



(ST) Z 



.8Tj„T B (t) 



f (ST) = P e 

If we now follow a procedure similar to that sketched in Ref. [T3|, we may write the relation 

B le (e B )=e-^f eB (e-( h ^f) e *) f , 



'ID 



where (h) f is the average of h(5T) over the distribution / f£T). ()10J) . By expanding the exponential 
within the angular bracket up to second order in its argument and after evaluating the average we 
obtain 



Bi e (e B ) ~ e 



-(h) f e B 



1 - -S^e 



2„2 



;i2) 



where we have introduced the abbreviation S 2 = 



It should be stressed that this result was 



obtained from a nonequilibrium thermodynamic analysis and expresses the effective Maxwell-Boltzmann 
factor in an universal form which only depends on the second moment, (5T 2 )f, of the temperature 
fluctuations. This corresponds to the low-E universal behavior discussed in Ref. [13 , where this 
universality allows for the definition of a universal parameter for any superstatistics and not only for 
Tsallis statistics. 



3 Hydrodynamics with fluctuating transport coefficients 

As mentioned in the introduction, fluctuations in intensive parameters such as the temperature, can 
not be described by the usual Onsager-Machlup theory, since it is devised to describe fluctuations in 
the extensive parameters only ;7j. For this reason, here we will describe the dynamics of a Brownian 
particle in the presence of large local temperature fluctuations ST(t) in its (r, u)-space, where r and u 
stand for the position and the velocity of a Brownian particle. More precisely, the state of the particle is 
defined not only by r and u, but also by the value of the parameter ST(t) entering in the single particle 
nonequilibrium probability distribution in the form P(f,u;t,ST). To accomplish the description, note 
that at the mesoscale the nonequilibrium probability distribution P obeys the continuity equation 

-P(f, u; t, 5T) + V- [uP] = -— ■ (Pv u ) , (13) 
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where V = J^, ^(r, w; t, 5T) is the corresponding streaming velocity in u -space. In this case, the 
explicit form of Eq. (pHIj) can be constructed by using Eq. (JJJ) and by using the following functional for 
the local nonequilibrium Kullback-like entropy s(f,t) [13 US EH 1201 CHI M 

s(r,t) = -k B [ P(r,u;t,8T)hi^- du, (14) 

J "le 

where P/ e characterizes a (local) equilibrium reference state, and is given by Eq. (J7J) for 5p = 0. It is 
important to stress that the entropy functional (|T4*|) will be used here as an irreversibility criterion that 
reduces to the Gibbs entropy in the case of equilibrium [TBI EH HJJ ■ 

By taking the time derivative of Eq. ()14|) and by substitution of the result into Eq. (|13|). one may 
calculate the entropy production a(f, t) of the Brownian gas through its natural evolution in time. After 
integrating by parts and by assuming that the currents vanish at the boundaries, one arrives at 

ds f d P 

— + V- J s = -k B / Pvf — In— du, (15) 
at J du Pi e 

where we have defined the entropy flow by 

J s = -k B J Pu [\n ^- - lj du. (16) 

For convenience we have assumed that the average temperature of the bath is constant. The right hand 
side of Eq. (|T5|) is the entropy production in the space of variables determining the state of the gas in 
the mesoscale, [JJ5], [2U| . 

Given the entropy production in ()15j) . one may impose the condition of a (locally) positive entropy 
production in the (f , ?i)-space and by following a scheme similar to that of nonequilibrium thermody- 
namics, in first approximation formulate linear laws between the current Pvn and its conjugated force 
J|ln-^-, [Uj, ^HJ- After taking into account Eq. (0), the linear law takes the form 

Pv a = (3h(5T)uP + -^fc (17) 

m ou 

where j3 is the corresponding coupling coefficient which has been assumed to be a scalar function for 
simplicity. Near equilibrium, this coefficient must reduce to the corresponding Onsager coefficient. Now, 
by substituting Eq. (fTTjl into f(13j) . we finally arrive at the following Fokker-Planck equation 

f) r) ( k T f)P\ 

-P (r>1M T) + V.[ S P] = ^. ( W maP + ^_), («) 

which describes the evolution in time of the nonequilibrium probability distribution of the Brownian 
particle. It is important to point out that the factor h(ST) breaks the equilibrium form of the fluctuation- 
dissipation theorem (FDT), thus implying that in the presence of large temperature fluctuations the 
usual form of the FDT is not valid. This result is due to the lack of thermal equilibrium between the 
system and the heat bath. 
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Once obtained the Fokker-Planck equation (|18j) it is now possible to derive an hydrodynamic-like 
hierarchy of equations for the moments of the probability distribution function P(f, u\ t, ST). As shown 
below, these equations constitute a fluctuating hydrodynamic description since they contain fluctuating 
transport coefficients. The evolution equation for the mass density field p(f; t, ST) = m J Pdu is simply 
obtained by directly integrating Eq. (JT5|l over the velocity u, yielding 

J^(r; t, ST) = -V ■ [pv(r; t, ST)} . (19) 

In a similar form, the evolution equation for the momentum field pv(r; t, ST) = m J uPdu can be shown 
to be 

p^-vif- 1, ST) + V-P k = -p/3* ■ v, (20) 
at 

where (3* = (3h (ST(t)) is a fluctuating friction coefficient. Here we have introduced the kinetic part of 
the pressure tensor P fe 

P^um^jiu-^-^, (2!) 



which obeys the following relaxation equation 



— P k + 2 
dt 



(3*f+Vv + \{V-v)i 



2k R T R i 

B -Ppl, (22) 



obtained by taking the time derivative of Eq. ()21|) and using (fTHj) . Here 1 is the unit tensor, the 
superscript s denotes the symmetric part of a tensor and we have neglected the contribution of the 
third and higher order moments (see Ref. [2H1 )• Eqs. (jl9|) - (}2*2*j) constitute a generalized fluctuating 
hydrodynamic equations of the Brownian phase when large temperature fluctuations drive the system 
out of equilibrium, breaking its thermal equilibrium with the bath j22j. It should be remarked that if 
the usual internal thermal fluctuations associated with the term proportional to ksTs were to be taken 
into account, a random heat and stress flows should be added to the elements of the pressure tensor 

P^ Then Eqs. (|2T)|) and ()22j) may be rewritten in the usual form of the Landau and Lifshitz fluctuating 
hydrodynamic equations, [22j. 



3.1 The limit of long times 

At sufficiently long times t/3* » 1, the description of the Brownian motion can be performed in terms 
of a diffusion equation for p(f; t, ST). This equation can be derived by neglecting the time derivatives in 
Eqs. (|2T)|) and (j2^|) . Then, Eq. (j2^|) leads to the following constitutive relation for the pressure tensor 

p*~h*EB ph -\6T% (23) 
m 

where, for simplicity, the Brownian contribution to the viscosity has been neglected j22l E] . It should be 
noted in passing, that precisely the same first two terms on the right hand side of Eq. (J18)) . which were 
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responsible for the breaking of the fluctuation-dissipation theorem in its usual form, are also responsible 
for the appearance of the correction h~ 1 (8T) on the right hand side of Eq. (J2HJ). Substitution of (|2*3|) 
into Eq.()20p leads to the following expression for the diffusion current pv 



pv 



D fluc (ST)Vp, 



(24) 



where the fluctuating diffusion coefficient has been defined as 



D fluc (ST) 



m/3 



h-\5T). 



(25) 



From Eqs. (|2*Hj) and (|24j). one may expect that when the fluctuations 8T are spatially non-homogeneous, 
then a fluctuating heat flow contributes to the diffusion flow. These results imply that in the long time 
limit under consideration and when large temperature fluctuations occur in the system, the diffusion 
coefficient of the Brownian particle is converted into a fluctuating quantity The evolution equation for 
the mass density p(r; t, 5T) in this limit is the obtained by substituting Eq. (|2*4*j) into (fH3j) yielding the 
fluctuating Smoluchowski equation 



which describes the evolution of the Brownian particle for long times. Eq. (|2fij) reinforces the fact that 
the fluctuation-dissipation theorem in its usual form is no longer valid. 

4 Regression laws in Brownian dynamics 

We can also derive the regression laws for both, the conditional averages of the velocity of the Brownian 
particles, and of the temperature fluctuations over a given volume element dV = dr. To accomplish this 
objective we must take the average of the corresponding Fokker-Planck equation that also incorporates 
the temperature fluctuation 5T , over df and then use the definitions of the conditional averages for 
the velocity and for the temperature 



d_ 

df 



■p(r;t,5T) = D fluc (6T)V 2 p, 



(26) 




(27) 



and 





The upper bar denotes a given initial state and 




(29) 
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Using this definition, from the corresponding Fokker-Planck equation a time evolution equation for 
P u (u, 6T; t) is readily obtained. If we insert the resulting equation into the time derivatives of Eqs. (J27j) 
and (|28|). after a partial integration similar to that of Sec. Ill, we arrive at the equations 

!^«> = -,S(HV( t ), (30) 

and 

^ = -^sr(t), (3D 

where e is an Onsager coefficient and we have assumed that (3 and e are independent of position. The 
nonlinear generalized force X ST (t) has been defined as 



and the modified friction coefficient (3 through the relations 

/?(5T°)if(*) = y uf3h{5T)P u dud{5T) 

= J uij ph(5T)P u d(5T)\du = j3(Ff) j uf(u,t)du, (33) 

where f(u,t) = f P u (u,ST;t)d(5T) is a reduced distribution function. By taking into account the ex- 
plicit expression for h(5T) and the fact that the average of 6T vanishes, one obtains [} — (3 (l — 7jk-(5T°) 2 
Note that the regression law for the temperature fluctuations is a non-linear one, and that it is coupled 
with that for u (t). 

From Eqs. (|3U|) and (|31|). it follows that mean regression laws are obeyed by the time-dependent 
value of the fluctuations occurring in the system after being averaged over an elemental volume. It is 
interesting to note that in the present case, when external factors induce temperature fluctuations, the 
average coefficients characterizing the dissipation in the system become modified by a term proportional 
to the square of the amplitude of the fluctuation. Eqs. (j3U|l and (}3~T]) may be interpreted as the gener- 
alization of the Onsager regression hypothesis in which the presence of large temperature fluctuations 
introduce corrections to the dissipation coefficients. 



5 DISCUSSION 

By using an irreversible thermodynamic approach we have derived an effective Maxwell-Boltzmann 
factor (EMBF) for a gas of Brownian particles, when there occur large fluctuations in the temperature 
of the heat bath in which they are suspended. The derived EMBF has the same form as the one the 
proposed in the literature [13j. Our derivation was carried out by using irreversible thermodynamic 
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arguments of the type used to evaluate the minimum amount of work that a system embedded into 
a heat bath needs to perform in order to modify the state of an external body [T3|. The EMBF has 
been expressed in a form that relates it to the superstatistics discussed in Ref. |T3]. The general 
point of view we have presented here indicates that this EMBF had its origin in a contraction of the 
mesoscopic space of variables necessary to completely describe the dynamics of the system. In this case, 
the superstatistics present in this mesoscopic stochastic variables have been eliminated by averaging 
over them. As a consequence of this contraction the extensive character of the description is lost. 

A hydrodynamic description was derived from a Fokker-Planck equation for the nonequilibrium 
distribution function of the Brownian system. The last equation was obtained by using a Kullback-like 
entropy functional which incorporates the complete distribution function as a reference state, and 
then averaging over the fluctuating velocities of the Brownian particle. The resulting description can be 
considered as a generalized fluctuating hydrodynamics scheme for the Brownian gas, that incorporates 
fluctuations of the intensive variables through random transport coefficients. Finally, by studying the 
appropriate space averages, we have derived non-linear regression equations for the time dependent 
average of the variables of the Brownian system. 

The results we have obtained in this paper, show that the description of far from equilibrium 
systems given in (r, u, <5T)-space lead, after contraction, to a non-extensive description. They also lead 
to a generalized hydrodynamics with fluctuating transport coefficients. Actually, since the transport 
coefficients depend on 8T, the linear relation between fluxes and forces in the (r, <5T)-space is in general 
lost, implying again a nonextensive description 
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